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Abstract 

We consider three two-level àtoms inside a one-dimensional cavity, interacting with 
the electromagnètic field in the rotating wave approximation (RWA), commonly 
used in the atom-radiation interaction. One of the three àtoms is initially excited, 
and the other two are in their ground state. We numerically calculate the propa- 
gation of the field spontaneously emitted by the excited atom and scattered by the 
second atom, as well as the excitation probability of the second and third atom. The 
results obtained are analyzed from the point of view of relativistic causality in the 
atom-field interaction. We show that, when the RWA is used, relativistic causality is 
obtained only if the integrations over the field freqüències are extended to — oo; on 
the contrary, noncausal tails remain even if the number of field modes is increased. 
This clearly shows the limit of the RWA in dealing with subtle problems such as 
relativistic causality in the atom-field interaction. 
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1 Introduction 



The propagation of electromagnètic signals has been a subject of investigation 
since the beginning of the quantum theory of the electromagnètic field. In 
the last years this subject has received much attention in the framework of 
rigorous proofs of relativistic causality in the atom-field interaction [1,2,3,4], 
of the recently observed slow light propagation [5] and of the possibility of 
superluminal propagation of light [6] . 
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Recently, the propagation and the scattering of a photon spontaneously emit- 
ted by an atom in a one-dimensional cavity has been investigated using the 
rotating wave approximation (RWA) [7]. This approximation consists in ne- 
glecting non-energy conserving terms in the Hamiltoniam, and the potential 
dangers of this approximation in terms of the causal behaviour of the physi- 
cal system considered are well known [8,3]. Quite frequently, the effect of the 
RWA seems to be compensated by the extension to — oo of the integrations in 
the frequency of the field modes, as in the original Fermi model description of 
causality in the excitation transfer between two àtoms, one of which initially 
excited [9,10,11]. Yet, this is in any case conceptually unsatisfactory for a sub- 
tle problem such as relativistic causality, and much effort has been dedicated 
to the inclusion of the counterrotating terms in order to have a rigorous proof 
of relativistic causality in quantum electrodynamics [4,12]. 

In this letter we consider the same system considered by Purdy et al. in [7]. 
This system consists of three two-level àtoms inside a one-dimensional cavity, 
interacting with the electromagnètic field in the RWA. We show numerically 
that the causal behaviour obtained in [7] (using the RWA) indeed derives 
from an extension to — oo of the frequency integrations; this extension, how- 
ever, has no physical justification since it makes the Hamiltonian unbounded 
from below, and this is in general a fundamental point for the causality prob- 
lem, as pointed out by Hegerfeldt [13]. We show that even if the Hamiltonian 
is bounded from below, noncausal terms in the field propagation are present 
when the rotating wave approximation is used; these terms do not vanish when 
the number of field modes is increased. This explicitly shows that the RWA 
is not an appropriate approximation for dealing with problems of relativis- 
tic causality in matter-radiation systems, and that the counterrotating terms 
should be included. 

Our system is the same as in [7]: it consists of three two-level àtoms, named 
1, 2, 3, in a one-dimensional cavity. The cavity has length L, with two parallel 
plates at x = and x — L, and the positions of the àtoms 1, 2, 3 are L/4, L/2 
and 3L/4, respectively. In the Coulomb gauge and multipolar coupling scheme, 
within dipole approximation and using the rotating wave approximation, the 
interaction of the three àtoms with the radiation field is described by the 
following Hamiltonian (with units such that H — 1) 

3 _ 3 

H = 12 u i S i + 12 u n a\a n + 1212 (gjna n S+ + g* n a ] n Sj) (1) 

j=l n n j=l 



where index n — 1,2,3,... denotes the field modes of the cavity and j = 1,2,3 
denotes the àtoms. a n , are the annihilation and creation operators of the 
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n-th mode, Sf, are the pseudospin operators of the atom j, and 



g jn = üj sin — — ; Qj = ^-j- (2) 



where Uj is the transition frequency of atom j and fij is its elèctric dipole 
moment. In the expression of the coupling constant the near-resonance ap- 
proximation u n ~ ojj has been used. 

We wish to stress that, because we are describing our system in the multi- 
polar coupling scheme, the field operator that we calculate, the momentum 
conjugate to the vector potential, is the transverse displacement field Dj_(r). 
This operator, outside the sources, coincides with the total elèctric field E(r). 
This is an essential point, because E(r) satisfies a retarded wave equation, and 
thus it is expected to manifest causal propagation in space. On the contrary, 
when the minimal coupling scheme is used, the field operator conjugate to the 
vector potential is the transverse elèctric field Ej_(r), which is not a retarded 
operator; in fact, the source term in the corresponding Maxwell equation is the 
transverse current density, which is not localized in space for atòmic systems 
[14]. 

Our initial state is the state with atom 1 excited, àtoms 2 and 3 in their 
ground states and the field in the vacuum state, that is | e,g,g,0k)- In the 
RWA the only states participating to the evolution of the initial state are, 
with an obvious meaning of the symbols, | g, e, g, n ), | g, g, e, 0„), | g, g, g, l n ). 
Thus we can write the general state at time t as the following superposition 



iP(t)) = d(t) | e, g, g, 0„) + c 2 (í) | g, e, g, n ) 

+ c 3 (í) | g, 9, e, 0„) + J2 b n(t) I 9, 9, 9, W (3) 



The Schròdinger equation gives the following set of coupled differential equa- 
tions for the coefficients 



àj(t) = -i yco jCj (t) + 9jnb n (t)j (4) 

bn(t) = -iL n b n (t)+j29%C j (t)) (5) 



We assume that the àtoms 1 and 3 have the same transition frequency, u>i = u 3 , 
and we indicate with 5 = oj\ — uji the detuning of atom 2 compared with àtoms 
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1 and 3. We also put L/c = 1, which means using as the unit time the time 
taken by the light to cross the cavity. 

Approximate analytical solutions of these equations have been obtained [15], 
using a method based on Laplace transforms [16]. We have obtained numeri- 
cal solutions of these equations. We integrate numerically the set of differen- 
tially equations (4,5) with the Adams-Moulton-Bashfort method [17,18]. This 
multistep method is an algorithm more sophisticated than the Runge-Kutta 
method, typically used for this kind of problems, allowing to obtain more ac- 
curate results and to use a much larger number of cavity modes. In order 
to facilitate comparison of our numerical results with the analytical results 
obtained in [15], we shall adopt the same numerical vàlues of the detuning 5 
and of the atòmic decay rates jj —\ Qj | 2 . Therefore, we use 71 = 1, 72 = 16, 
73 = 256 and 5 = 4 (in our units). 

We calculate the expectation value of the square of the elèctric field, which 
is proportional to the elèctric part of the field energy density, given by (zero- 
point terms have been neglected) 



Our first numerical calculation involves a set of equally spaced field modes, 
symmetric with respect to the resonance frequency of the first atom. This 
simulates an integration extended from —00 to +00, therefore a field Hamilto- 
nian not bounded from below (as that used in [7]). In this paper we report the 
results of two numerical calculations: the first uses 10 4 modes, 5000 above and 
5000 below the atòmic frequency; the second uses 2 • 10 4 modes, 10 4 above the 
atòmic frequency and 10 4 below. We then compare the results obtained with 
those of an analogous numerical calculation in which the field modes are not 
symmetric with respect to the atòmic frequency, and only the upper cut-off is 
increased with increasing the number of the field modes; this second case sim- 
ulates a field Hamiltonian bounded from below. This comparison allows us to 
understand the role of the RWA in the behaviour of our system, in particular 
from the point of view of relativistic causality. 

Fig. 1 shows {E 2 (x, t)) at t — 0.25 when 10 4 modes are used with a symmetric 
distribution around the transition frequency. A front in the propagation of the 
energy density is evident at x — 0.5, as expected (the atom that emits the ra- 
diation is located at x — 0.25); however, a zoom of Fig. 1 in the neighbourhood 
of x = 0.5, given in Fig. 2, shows the presence of tails for x > 0.5, blurring 
the front. Fig. 2 also shows the same expansion around x = 0.5 when a larger 
number of modes is used (2-10 4 ), again with a symmetric configuration around 
the atòmic frequency. An improvement of the behaviour from the point of view 
of causality is evident, with a manifest decrease of the tails for x > 0.5. We 




(6) 
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have also obtained similar results for different times, showing that the front 
gets sharper when the number of field modes is increased, if the freqüències of 
the field modes are symmetric around the transition frequency of the atom. 
This suggests that the expected causal behaviour is indeed approached in the 
limit of an infinite number of field modes symmetrically distributed around 

Lüi. 

The behaviour is quite different when the number of field modes is increased in 
such a way that only the upper cut-off frequency increases but the lower cut- 
off frequency remains fixed (in this case the freqüències of the field modes are 
not symmetric with respect to the atòmic transition frequency); this situation 
is intended to mimic the case in which the field modes extend from to oo, 
and no extension to — oo is performed. Fig. 3 shows a zoom around x = 0.5 of 
the square of the elèctric field at time t = 0.25 for 10 4 , 2- 10 4 and 3-10 4 modes. 
The figures indicate that, as the number of modes is increased, the envelop of 
the oscillations does not approximate a sharp front. The remaining tail gives 
a noncausal behaviour. This makes quite evident that the use of the rotating 
wave approximation with an Hamiltonian bounded from below does not give 
causality in the propagation of the electromagnètic fields, because noncausal 
tails persist, consistently with Hegerfeldt's theorem. 

Similar conclusions are reached by calculating | cz(t) | 2 , that is the excitation 
probability of atom 3. Causality requires that this probability should vanish 
for t < 0.5. Fig. 4 shows our numerical results, and Fig. 5 a zoom around the 
causality time t = 0.5 (with 10 4 and 2 • 10 4 modes), both with a symmetric 
configuration of the mode freqüències. We note that a sharper front is obtained 
when the number of the modes is increased. However, Fig. 6 shows the result 
for the non-symmetric configuration, in which only the upper cut-off frequency 
is increased with the field modes. It is evident from Fig. 6 that in this case 
the noncausal tails for t < 0.5 do not decrease on the average as the number 
of the modes is increased. 

We wish to conclude by stressing that our results clearly show the reason why 
recent results in the literature have obtained a causal behaviour of the atom- 
field interaction in a cavity within the rotating wave approximation. The point 
is that at the same time the frequency integration over the field modes was 
extended to — oo which permits to escape the conditions set by Hegerfeldt's 
theorem by making the Hamiltonian unbounded from below, but which is 
physically unacceptable. In this paper, we have considered three two-level 
àtoms, one excited and two in the ground state, inside a one-dimensional cav- 
ity, interacting with the electromagnètic radiation field in the RWA. We have 
calculated numerically the energy density of the elèctric field spontaneously 
emitted by the excited atom and scattered by the second atom, as well as the 
probability of excitation of the second and third atom. We have shown that, 
without the (arbitrary) extension of the field freqüències to — oo (frequently 
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used in the literature), and which has no physical basis, noncausal tails are 
present both in the field propagation inside the cavity and in the atòmic ex- 
citation probabilities, even when the number of the modes is increased. This 
underlines the potential dangers of the rotating wave approximation. In a 
forthcoming paper, we will explicitly show that the correct inclusion of the 
counter-rotating terms of the Hamiltonian, allows to obtain a better causal 
behaviour without any need of extending the frequency of the modes to — oo. 
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Fig. 1. Average value of the square of the elèctric field E 2 (x,t) in arbitrary units 
at t = 0.25, with 10 4 field modes symmetrically distributed around the atòmic 
transition frequency. 
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Fig. 2. Average value of E 2 (x, t) in arbitrary units at t = 0.25 in the neighbourhood 
of x = 0.5. The dashed line is with 10 4 modes of the field and the continuous line 
with 2 • 10 4 modes (both with a symmetric mode distribution around the atòmic 
transition frequency). 
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Fig. 3. Average value of E 2 (x, t) in arbitrary units in the neighbourhood of x = 0.5 
at time t = 0.25, with 10 4 modes (dashed line), 2 • 10 4 (dotted line) and 3 • 10 4 
(continuous line). The distribution of the field modes used here is not symmetric 
respect to the atòmic transition frequency. 
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Fig. 4. Excitation probability of atom 3, with a configuration of the field modes 
symmetric around u\. 



11 




Fig. 5. Zoom of Fig. 4: excitation probability of atom 3 with a symmetric configu- 
ration of the field modes, around the causality time t = 0.5. 
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. 6. Zoom of the excitation probability of atom 3 around the causality time 
0.5, using a non-symmetric distribution of the field modes. 
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